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We study a generic model for the polarisation and motility of self-propelled soft objects, biological 
cells or biomimetic systems, interacting with a viscous substrate. The active forces generated by 
the cell on the substrate are modelled by means of oscillating force multipoles at the cell-substrate 
interface. Symmetry breaking and cell polarisation for a range of cell sizes naturally “emerge” 
from long range mechanical interactions between oscillating units, mediated both by the intracel¬ 
lular medium and the substrate. However, the harnessing of cell polarisation for motility requires 
substrate-mediated interactions. Motility can be optimised by adapting the oscillation frequency to 
the relaxation time of the system or when the substrate and cell viscosities match. Cellular noise 
can destroy mechanical coordination between force-generating elements within the cell, resulting 
in sudden changes of polarisation. The persistence of the cell’s motion is found to depend on the 
cell size and the substrate viscosity. Within such a model, chemotactic guidance of cell motion 
is obtained by directionally modulating the persistence of motion, rather than by modulating the 
instantaneous cell velocity, in a way that resembles the run and tumble chemotaxis of bacteria. 


I. INTRODUCTION 

Cell motility on solid substrates (crawling) and in fluid 
environments (swimming) are usually regarded as being 
based on different physical principles and studied inde¬ 
pendently. Swimming often relies on the beating or ro¬ 
tation of protrusive appendages (flagella or cilia) Bi, 
while crawling generally relies on the protrusive and con¬ 
tractile forces generated by the acto-myosin cytoskeleton, 
and on cell-substrate adhesion [I[. Remarkably, several 
cell types that undergo major shape changes while crawl¬ 
ing, such as amoebae and neutrophils, are also able to 
swim in Newtonian fluids ii- Furthermore, crawling 
cells are sensitive to the stiffness (elastic response) [B- 
0, but also to the viscosity [1, @ of the substrate on 
which they are crawling. This suggests that crawling 
and swimming share common underlying physical prin¬ 
ciples and that insight on eukaryotic cell motility may 
be gained by studying self-propelling soft objects in fluid 

environments [ 13 - 

Cell motility requires polarisation, namely the break¬ 
ing of front-back symmetry of the cell. This can be trig¬ 
gered by external gradients of either biochemical (chemo¬ 
taxis) [ll| or mechanical {e.g durotaxis) [13 nature, or 
it can occur spontaneously. Spontaneous cell polarisa¬ 
tion is of particular interest to understand the spatio- 
temporal correlations inside and outside the cell. Active 
gel theories [Tsl - fl^ have shed some light on the physics 
of spontaneous cell polarisation, but much remains to be 
learnt about how coordination is achieved at the scale of 
the whole cell, and how coordinated motions are affected 
by the properties of the extra-cellular environment. 

Here, we devise a generic theoretical framework which 
enables us to address in a unihed manner three different 
fundamental aspects: polarisation, motility and chemo¬ 
taxis. Our model soft cell (henceforth simply ‘cell’), 
which aims at describing active droplets [ISl - Hh . H^ , bio¬ 
mimetic self-propelled systems and cells, is a viscous or 


an elastic body interacting with a substrate by means 
of localised forces distributed over the cell-substrate con¬ 
tact area. We restrict ourselves to the case where the 
substrate is a Newtonian fluid, so that force transmis¬ 
sion between the cell and the substrate can be achieved 
by enforcing a no-slip boundary condition. Our findings 
also provide insight into the polarisation of cells adher¬ 
ing to elastic substrates, although crawling in that case 
requires additional assumptions concerning the dynamic 
of creation and destruction of adhesion sites. 

The cell cytoskeleton commonly displays quasi- 
periodic spatio-temporal patterns and fluctuations of ac¬ 
tivity [20M^. These oscillations are observed in both 
adherent and crawling cells l26l-[^ even in the 

absence of specific adhesion receptors . Motivated by 
that, we model the distribution of cellular forces as os¬ 
cillating force-multipoles at the cell-substrate interface. 
Spatio-temporal correlations among the oscillating forces 
may spontaneously emerge from mechanical interactions, 
affecting cell polarisation and cell motility. To show this, 
we first study the motility of a cell exerting forces on 
the substrate with prescribed time-dependence. The re¬ 
sulting cell speed has non-monotonic trend as function of 
the substrate viscosity or of the frequency of the oscillat¬ 
ing units. We then study spontaneous symmetry break¬ 
ing associated to cell polarisation and the transition from 
non-motile to motile states as a dynamical process emerg¬ 
ing from the phase-locking of the oscillating elements. 
This process occurs only for a range of cell sizes and is 
permitted by long-range interactions mediated by both 
the substrate and the cytoplasm. To investigate chemo¬ 
taxis, we introduce a minimal coupling between chemo¬ 
tactic gradients and the force distribution. Finally, we 
discuss the persistence of cell motion in the presence of 
noise (e.g. cellular noise) and we obtain testable predic¬ 
tions regarding the substrate-dependent statistical prop¬ 
erties of the cell trajectories. 
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FIG. 1. (Color online) A motile cell exerting traction forces 
on a flat viscous substrate is modelled as two coupled fluids: 
a layer of thickness fi and viscosity rjc (the “cell” on top, in 
green) lying over an infinite bulk fluid substrate of viscosity 
rjs (at the bottom, in light blue). Cytoskeletal elements gen¬ 
erating traction are modelled as oscillating force multipoles 
(two dipoles in the sketch) distributed along a line. 


semi-infinite fluids, f) —> oo. The interfacial velocity field 
at a distance r from a unique interfacial point-force f is 
“ 2-K{iu+r^a) r^' effect of several point-forces is 
additive thanks to the linearity of Eq. O- 

To study the emergence of cell polarisation we consider 
time-dependent and periodic force distributions consist¬ 
ing of two identical units A and B, labelled with index a 
in the following, that on average over a period are mirror 
symmetric with respect to the cell centre (FiglT|). For 
simplicity, all forces are distributed along a line. Each 
unit is made of particles with coordinates x“ = 
subjected to forces F“ = F^x directed along the x axis. 
Each unit centre is at Cq, = 2/fV^„x“, the cell center 
at c = (c^ -b Cb)/ 2, and the separation between units is 
r := Cb - ca- 

III. SUBSTRATE VISCOSITY AFFECTS THE 
CELL SPEED 


II. MODEL 

A. Derivation of the cell speed 


To begin with, we model the cell as an incompressible 
fluid layer of thickness f) and viscosity rjc lying over an 
incompressible bulk fluid of viscosity r]s (Fig[T|) at low 
Reynolds number. The fluids are immiscible and satisfy 
the Stokes equation [s^: 

= 0; V • = 0, (0 < z < ()) (1) 

= 0; V • = 0, (-oo < z < 0) 

where and are the velocity and pres¬ 

sure fields in the fluid layer and the substrate, and z the 
direction normal to the interface. The fields and 
are coupled by means of conditions at the interface z = 0. 
In the following, we assume no-slip boundary conditions 
at the interface, and neglect normal 

displacements of the interface, = 0. 

Similar results are obtained if no-slip is enforced only at 
the location of the active force generators, see Appendix 

1X1 

Protein complexes responsible for the traction forces 
at the cell-substrate interface are modelled as N disks of 
finite radius a and negligible height e —0 distributed at 
positions r„ (n = 1,.. . N) at the interface between the 
two fluids via F{r, z = 0) = r'n<5(r~rn). By virtue 

of Newton’s third law, the cell exerts no net force on the 
substrate, X]n=i = 0, and the local stress balance at 
the interface reads 

where denotes the in-plane component of the veloc¬ 
ity for c = s, c, see Appendix lA II The solution of 
Eqs.ffll can be obtained using Fourier transforms, see 
e.g. [^,1^, as explained in Appendix I A II We focus 
here on the analytically tractable limit of two coupled 


Interactions between the point forces may lead to the 
motion of the force distribution with respect to both the 
cytoplasm and the substrate. Only the interactions me¬ 
diated by the substrate can lead to displacement of the 
whole cell. Since the cell’s radial boundary is not in¬ 
troduced explicitly in our model, the cell net velocity is 
obtained by subtracting the cytosol motion to the mo¬ 
tion of the force centre c. This leads to a force balance 
Cc = J2a 'lln explained in detail in Appendix IB II 

Here, 5“ = Csv(x“) represents the substrate-mediated 
traction force, and v(x“) := v('’’)(x“) = v(®)(x“) is the 
velocity at the location of disk (a, n) due to the motion 
of all remaining disks. The total and substrate-related 
drag coefficients are C, and Cg. In the limit t) —>■ oo for a 
disk of radius a : (s = and ( = ^[rjs + ric)a [s^ . 

We first consider the case where each unit is an os¬ 
cillating force dipole made of two particles: n = 1,2 
with one force scale = — Ef := fa and one length 
scale La = xf — xf- With no loss of generality we write 
La = la + da where la is a constant and da describes 
time-dependent deformations, see FiglT] We parametrize 
the amplitude and force of the oscillating dipoles as: 

da = RaCOs{ujt + 4>a); fa = -gaSin{ujt + 4>a)- ( 3 ) 


The average cell velocity is obtained from the net dis¬ 
placement of the cell centre over an oscillation period. 
For the pair of identical dipoles {{g,R,l}A = {gJ^Rfs) 
shown in Fig[T] for f) —>■ oo and to lowest order in la/r, 
it has a simple expression derived in Appendix IB 41 


=(V’) := 


T 


dtcif) 


gRP^ 


iVs + VcY 


sin Ip (4) 


where ip := (pB ~ (pA is the phase difference between the 
oscillators. Geometrical details of the force distribution 
are contained in S = 3f;f[4-|- {R/lY cos'?/;], see Appendix 
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IB 4dl Replacing the viscous cell with an elastic cell of 
elastic modulus yields a similar result, with (*w) in 
place of r]c, see Appendix lA 31 Appendix IB 21 Appendix 
|B4^ 

The relationship between forces and displacements 
is obtained from the force-balance equation: Cdo, = 
2fa +Ta- Here = ^[v(x“) — v(x 2 )] represent non¬ 
local interactions, propagated both by intracellular and 
extracellular media. For small amplitude oscillations: 
Ra la < r, interactions may be neglected to lowest 
order and the force-amplitude relation is Qa = {C^Ra)/‘2- 

The active forces (typically actin polymerisation and 
actomyosin contraction for crawling cells) may impose 
either the force or the displacement scale. Eq.d!]) shows 
that the migration speed is qualitatively different for im¬ 
posed displacement {vc ^ rjs/{rjs +ric) x au}{Rl)'^/r^) and 
imposed force {vc ~ Vs/iVs + Vc)^ ^ f^e 

latter case, the velocity presents a maximum when the 
substrate and cell viscosity are similar (see Figl^Ka) - the 
value of the optimal ratio depends on "0). This biphasic 
behaviour shows an interesting analogy with the biphasic 
velocity of cells crawliM on elastic substrates in response 
to substrate stiffness p-Q- The speed of an elastic cell 
can be optimised by tuning the frequency of the oscil¬ 
lating force units close to the inverse relaxation time of 
the system: a; ~ 2^, see Figl^Kb). This suggests that 
cells may adjust their oscillations rate to the mechani¬ 
cal properties of the environment for optimal motility. 
By symmetry, a fluid cell with oscillating disks perma¬ 
nently bound to an elastic substrate experiences a net 
cytoplasmic flow. However, harnessing this flow for cell 
migration on an elastic substrate requires additional hy¬ 
pothesis concerning disk attachment to and detachment 
from the substrate, or disk creation at the cell front and 
destruction at the rear. 

Using typical values of parameters obtained from ex¬ 
periments we estimate Vc ~ I0{fj,m)/hr. This is smaller 
than typical cellular speed, a common drawback for 
swimmers, made up of point-like particles, undergoing 
small amplitude strokes [l^- Higher speeds can be ob¬ 
tained by increasing a/r and la/r beyond the validity of 
the analytical results, Eq.(|4]). This expression is however 
very valuable, as it shows that net migration requires 
that the oscillators phase lock at U ^ 0(7r), like inter¬ 
acting dumb-bells in bulk fluids [s^ This breaks 

time-reversal symmetry, a requirement for motility at low 
Reynolds number (“scallop theorem” [3^1. 


B. Estimated value for the average migration speed 



FIG. 2. (Color online) (a) Average migration speed (Eq.(|4])) 
as function of the substrate viscosity ris. The force distribu¬ 
tion is shown in Fig[T] with ip = 7r/2, for either fixed force 
(continues black line) or fixed displacement (dashed red line). 
The former is measured in units of (9g^l^)/{6‘iTTri^ujar^) and 
the latter in units of FI?)/{'kt'*'). rjs is measured in units 
of r^c. (b) Average migration speed variation with the driving 
frequency for an elastic cell of elastic modulus pc- Velocities 
for fixed force (dashed dotted black line) and fixed displace¬ 
ment (dotted red line) are measured respectively in units of 
{9g^?)/{64'KrjsiJ,cCLr‘^) and in units of {apcR'^?)/(gsTTr^) while 
u is measured in units of Hc/Vs, the viscoelastic frequency. 


cellular measurements give ~ lOOPa at i/q ^ IHz, 
see [ 3 ^. The associated cell viscosity is G^/w = ^Pa s 
(a highly viscous regime: for comparison water at room 
temperature has a viscosity ~ 10“^Pa s - four orders 
of magnitude below). For simplicity we take a substrate 
with equivalent viscosity, i.e. we consider conditions of 
‘viscosity matching’ between the cell and the substrate. 

We next address the relation among force and defor¬ 
mation for the oscillator. We consider forces g ^ 0.5nN 
that are varying with frequency w ^ 2'iT{rad)/s drag¬ 
ging a disk of size 0.5fim. The drag coefficient is C = 
[16a(?7s -I- ? 7 c )]/3 ~ (10aG^)/a;. The corresponding oscil- 
_ 2^ 
uiC 


lation amplitude is P = ^ ~ loxo^s^mmoPa ^ 2/rTO. 


The average migration speed can be calculated us¬ 
ing the formula valid for fixed deformation amplitudes, 
reported above. In addition to the previous values, 
we consider I ^ 2R ^ 4fim and r ^ 10/rm (as well 
as ‘viscosity matching’ 77 ^ = rjc) from which Vc 
f' ~ [lOpm)/hr. Under same conditions 

of lengths, and matching between cell-substrate moduli 
77 s ~ ^ , the elastic cell has the same speed of the viscous 

iz>2i2 3 

cell Vc ^ ^ 4 *^ 




IV. CELL POLARISATION, MOTILITY AND 
CHEMOTAXIS EMERGE FROM 
SYNCHRONISATION 


Probing the cell with periodic stress or strain of well 
defined frequency ui allows one to infer the viscosity of 
the cytosol rjc, from measured values of the shear loss 
modulus Gj, by means of the relation r/c ~ Cl^/w. 

To estimate values of migration speed in our model, 
we consider oscillation frequencies vq ^ IHz, (a lower 
bound for swimming cells [ 2 |) i.e. uj = 2'K(rad)/s. Intra- 


We now show that polarisation and migration of the 
cell may spontaneously emerge as the result of dynamical 
interactions and synchronisation of the force generators. 
The amplitudes Ra(t) and gait), and the phase (pait) 
are now slowly varying quantities : constant over an os¬ 
cillation period T = 2'k/uj, but varying over l ong er time 
scales as a result of intracellular interactions 1^ . A 
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simple model (40j leading to self-sustained oscillations is 
the evolution equation for the forces, 

fa = -^[—}Cada + A4/ct[l — tjd^] -|- Ad^] (5) 


combined with the force balance equation = 2^a+Ia 
(for a = B). These equations are equivalent to a van 

der Pol-Dufling oscillator E3, El that exhibits super¬ 
critical Hopf’s bifurcation in a wide range of the pa¬ 
rameter space. Similar models emerge from the collective 
dynamics of molecular motors, see e.g. [H, Here, 
Ka sets the oscillation frequency, M. determines the in¬ 
stability threshold and ct is a stabilising term. These 
parameters must be strictly positive for a stable limit 
cycle to exist. A 0 describes non-isochronous oscil¬ 
lations [i^, and determines whether the oscillation fre¬ 
quency increases {A < 0) or decreases {A > 0) with 
increasing amplitude Ra- This non-linear coupling may 
for instance result from the mechano-sensitive kinetics of 
sub-cellular constituents. 

In the absence of chemotactic bias, the parameters /C, 
a and A are identical at both ends A and B of the cell. 
We implement chemotaxis by writing K.^ = 1C + iCpa 
where pa = p{ca) is the density of chemoattractant at the 
centre Cq, of unit a. This description is consistent with 
recent studies on chemotaxis [^, where pseudopods at 
the cell edges display wave-like patterns and chemoat¬ 
tractant changes the rate of internal processes (here the 
oscillation rate, see below). 

The dynamics of a single oscillator is seen neglecting 
interactions X„. Combining Eqs. (I515|) with Qda = 2/a, 
and averaging over the (fast) period T yields Ra = 
A4i?a[l — cF{RaY 1^, showing that the amplitude sat¬ 
urates at a stable value 94 = and the phase 

rotates at slow frequency ^a = —JCa + .494^. To study 
synchronisation of distant oscillators, we assume that the 
interactions Ig induce small deviations from the limit 
cycle d^, see for details, yielding a phase equation: 


^ H-I-sgn[^] —(C-I-cos V’) sin')/). (6) 

To 


Here fl = ~ P^) the chemotactic bias, the 

relaxation time tq results from the long-range interac¬ 
tions between oscillators, and U contains geometrical de¬ 
tails of the force distribution. Eq.® resembles Adler’s 
equation 1^, previously used to describe synchronisa¬ 
tion dO, l46j| and swimming of algae having intrinsic 
front-back asymmetry [El - 13 iir low Reynolds num¬ 
ber fluids. The additional term sin-!/)cos/; is associ¬ 
ated to force multipoles. For two dipolar units (FiglT|); 
To = (7r7WCr3)/(2|A|944a) and U = 2P/^'^. 

It is instructive to write Eq.® as /; = 
in terms of an effective potential V with minima corre¬ 
sponding to (meta)stable phase-locking between the two 
oscillating units (FigO]). Motility requires ijj 
see Eq.dH). In the absence of chemotactic bias (H = 0), E 


presents a single minimum (modulo 27r) if |C/| > 1 (con¬ 
tinuous blue curve in FigI3J)). A bistable systems with 
■0 0(7'') spontaneously emerges from the long-range in¬ 

teractions if A > 0 and \U\ < I, with two symmetric sta¬ 
ble states corresponding to spontaneous cell polarisation 
and motion (dashed-dotted green curve in FiglSJo), sepa¬ 
rated by an effective energy barrier AEq = (1 ~ |t^l)^/2- 
Moderate gradients of chemoattractants (0 < ro|H| < 1) 
introduce a bias that displaces the single minimum (|17| > 
1 ) or favours one of the two stable states (|C/| < 1), (con¬ 
tinuous blue and dashed-dotted green curves in FigI3Jl), 
thus directing motility. 

The pair of dipolar oscillators considered so far (FigH]) 
has U > 1 and does not exhibit a stable state with bro¬ 
ken symmetry. Spontaneous cell polarisation is possible 
if the force distribution in each unit is itself polarised, 
the whole cell conserving mirror symmetry. In Fig]^, 
we show an example where each unit consists of two 
dipoles (/ and //), separated by a distance /. To keep 
the number of parameters minimal, dipole I is chosen to 
be a scaled version of dipole II: li/ln = 947 / 94 // = k 
with Ijj = = 94. The two dipoles oscillate at 

the same frequency uj and with a fixed phase differ¬ 
ence; oscillator / is in opposition of phase with oscilla¬ 
tor II, which satisfies Eq.([5]). This way, the two units 
are characterised by a single dynamic al p hase differ¬ 
ence 0 satisfying Eq. ®. As shown in [4d[ |. the migra¬ 
tion speed for this force distribution is still given by 
Eq.® with 5 = ^{4(«:^ — I)^ -|- + l]^cos0} 

and retains the qualitative trend shown in FiglJJ More 
interestingly, phase locking at 0 0= 0(7r) occurs when 

\U\ = \ < 1 (where 5 := i/r) and A > 0. 

Correspondingly, A = 


V. NOISY OSCILLATIONS RESULT IN 
DIRECTED OR PERSISTENT MOTIONS 

Uncorrelated noise, associated e.g. to the stochas¬ 
tic nature of actin polymerisation and motor pro¬ 
teins activity, influences intracellular synchronisation. 
Adding a random effective “force”, v, to Eq.dH): 0 = 
^E(0) -I- v\, results in fluctuations around a sta¬ 
ble minimum and occasional jumps between two minima 
of E leading to phase slips [H, H^. When \U\ > 1, 
jumps give a phase slip of 2/7 and do not change the 
cell velocity. Velocity fluctuations are controlled by lo¬ 
cal fluctuation of 0 around the minimum 0*, estimated 
to be of order ((0 — 0*)^) ~ A/9|E(0*) ^ A/||17| — 1|, 
where A is the intensity of the noise, with correlation 
{v{t)v{t')) = 2Aroi5(t — t'). The cell performs a ran¬ 
dom walk along the x axis in the absence of chemotactic 
gradients, or a random walk with drift with such gradi¬ 
ents, with a diffusion coefficient of order D ~ {v'^)to ~ 
{dpVc{ip*))'^ToA/{\U\ — 1) (continuous blues curves in 
Figl3t,e). When \U\ < 1, the two stables states 0± 
(with sin ip- = — sin 0_|_ ) have non-vanishing velocities 
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FIG. 3. (Color online) (a) Force-distribution corresponding 
to two polarised oscillators at either edge of the cell, a system 
that can show spontaneous symmetry breaking (|C/| < 1). (b- 
e) Synchronisation landscape and cell trajectories obtained 
from numerical simulations (see [4^ for details). All pan¬ 
els have A > 0 and either U > 1 (continuous blue curves) 
or \U\ < 1 (dashed-dotted green curves), (b) Effective cell 
polarisation landscape V{tp), and (c) trajectories of the cell 
centre c{t) with time in the absence of chemotactic gradient 
(fl = 0). In the presence of chemotactic gradient, (d), the 
effective landscape V{tp) is asymmetric and (e) the cell tra¬ 
jectories show a bias. Insets in (c,e) are trajectories in the 
absence of noise. In this case, the direction of motion for 
\U\ < 1 is arbitrary, and may be away from the chemotac¬ 
tic gradient. All parameters are the same for the continuous 
blue curves and dashed-dotted green curves except r (double 
for the dashed dotted green curves) resulting in [/ = —0.6 
(dashed-dotted green curves) and U = 2.5 (continuous blue 
curves). 


duced by chemotactic gradients is mostly due to a bias in 
the “run” time 7±, due to the asymmetry of the potential 
barrier {AV- > AV+ in FigI3}?), and only slightly to di¬ 
rectional bias in the cell velocities (|z;c(V’-)| > kc(l/'+)|)- 
This chemotactic strategy, reminiscent of the one of bac¬ 
teria [i^, has to some extent been observed for Dic- 
tyostelium (d^. [5ll|. 

The value of |17| is critical to spontaneous cell polari¬ 
sation, the persistence of cell motion, and the cell chemo¬ 
tactic response. Factors that affect the value of U will 
thus qualitatively affect the cell’s trajectories. For a given 
force distribution in each unit at the cell edges, increasing 
the cell size r monotonously decreases U. The condition 
|C/| < I corresponds to a range of cell size, and we thus 
predict that very small or very large cells (compared to 
the size of one force unit) should not show spontaneous 
polarisation. Another remarkable result is that in the 
fixed force regime (when 91 oc 1/C), the parameter f/ ~ C^ 
increases with increasing substrate viscosity. We further 
predict that not only is the instantaneous cell velocity 
affected by the substrate mechanics as depicted in Fig 121 
but also the persistence of cell motion, and the ability for 
spontaneous cell polarisation in the absence of external 
cues, could be strongly impaired under high substrate 
mechanical resistance. 
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Appendix A: Flow at the cell-substrate interface 
1. Two coupled fluid layers 


Vcifp-) = —Vc{'ip+) (Eq.Q). The cell’s trajectory is thus 
a succession of persistent motions separated by random 
changes of direction, resembling the run and tumble mo¬ 
tion of bacteria [^. The persistence time of directed 
motion is the first passage time between two consecutive 
jumps: 7± (X (^ . Typical trajectories without 

and with chemotactic gradients, obtained from numeri¬ 
cal simulations (see for details), are shown in Figl3J:,e 
(dashed-dotted green curves). 


VI. DISCUSSION 

Our model captures within a single framework two dis¬ 
tinct chemotactic strategies, namely the case where the 
gradient polarises an otherwise quiescent cell (|17| > I) 
and the case where the gradient biases a preexisting po¬ 
larised state showing transient directed motion (|f/| < I). 
Interestingly, in the latter case the directional bias intro¬ 


In this Appendix we outline the procedure to obtain 
the analytical solution of Eq. (HD considered in the main 
text. Related problems were studied previously also by 
other authors, see e.g. ref 

We consider two planar fluid interfaces, at position 
z = 0 and at position z = fj where £ is a unit vector indi¬ 
cating the direction normal to the interface. To proceed 
it is useful to decompose the velocity appearing in Eq. (HD 
v‘’(a;,2/,z) = Vy(x,?/, z) + vl{x,y, z)z in terms of the in¬ 
plane components v|(a;,?/,z) and orthogonal component 
vl{x, y, z) where the label <; = c, s indicates either cell or 
the substrate. A similar decomposition can be done for 
the gradient: V = Vy + where Vy = (^; ^) has 
components only in the (x, y) plane. Einally the Lapla- 
dan becomes = V| -b 

The boundary conditions at the interfaces are as fol¬ 
lows: no motions along z-axis, = 0, and 

continuity of the tangential component of the velocity 
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at the interface; no-stress condition at the cell-upper 
solvent interface (obtained neglecting solvent viscosity) 
? 7 c[V||ui'^^ -I- = 0. At the cell-substrate inter¬ 

face, z = 0, the conditions are no z-motions of the inter¬ 
face, • z\z^o = 0 and • z\z=o = 0; continuity of 
the tangential component, = V||®^ j^^o; stress bal¬ 
ance, -f = 

—F{r,z = 0). As explained in the main text, it suffices 
to focus on a single interfacial point force F{r, z = 0) = 
f(5(r). 

We search a solution using Fourier transforms. By ex¬ 
ploiting the geometry of the system, we decompose the 
Fourier vector k as k = q -|- pz where q lies in the plane 
{x, y) orthogonal to direction z. Thanks to this decompo¬ 
sition, a generic vector field E(a;, y, z) in three-dimensions 
in presence of the planar interface at z = 0 can be written 
as 

(Al) 


We proceed by taking the Fourier transform of the ve¬ 
locity components. The condition of no motion in the 
z-direction at the interfaces means (q, 0 ) = 0 and 
'Di‘^^(q, h) = 0 so the various stress conditions only in¬ 
volve terms V||‘’^ (q, z) for c = c, s. 

We now make an ansatz concerning the z-dependence 
of the Fourier components. For the flow in the sub¬ 
strate z < 0 we pose V||'*^(q, z) := S(q)e'?^. For the 
intracellular flow, in the region 0 < z < (}, we pose 
V||'^^(q, z) = C-(.(q)e“'J^ -I- C_(q)e'^^. The three terms 
C±, S can be determined using the boundary conditions 
for the stress and velocity given before. 

As we have three unknowns, it suffices to consider three 
of the four equations for the tangential velocities dis¬ 
cussed above. We choose the continuity of the tangential 
velocity at at z = 0, C+-|-C_ = S and the two equations 
for the interfacial stresses at z = 1 ): 0 = r]c[—C+e~'^^ + 
C_e'^'’]; at z = 0: gpc)—C+ + C_] — qr]s^ = —f- From 
the second equation we obtain C 4 . = Using this 

relation in the first equation we get S = C_(l -|- 
Inserting these relations in the third equation, we deter¬ 
mine C_,C+ and finally the flow at the cell-substrate 
interface z = 0 


v(*)(q,0) = S = C+-l-C_ 


f 

g{ 7 ?c[tanh(()q)] -|- 77s} 


The inverse Fourier transform of S(q) has two analyt¬ 
ically tractable limits ; h —i' 0, corresponding to a thin- 
film, and f) — 00 , representing two semi-infinite coupled 
fluids. We focus on the latter where [tanh([i( 7 )] = 1 and 


= v^|'’-'(r,0) = 


is), 


{2Tr)‘^ {r]s + Vc) 


(fq-e 


2Tr{r)s + rjc) Jo 


dqJoiqr) = 


27r(?7s -f ric)r 


(A2) 


Here we used f d^q^e *'>■'■ = J^dqqJ^^dO^e ^rcos9 _ 

2tt dqJo [qr) where Jo indicates the modified Bessel 
function of the first kind of zeroth order • Eq. (IA2I) is 
the expression reported in the main text and represents 
the disturbance in the flow field at the interface of two 
semi-infinite fluids, generated at position r by a point 
force placed at the origin, and r := ||r||. We can write 
this relation as v(r) = H{r)i where H{r) is the Greens’s 
function describing the flow at the interface. Thanks to 
the linearity of the equations describing the fluids, the 
effect of N point forces is obtained by superposing the 
single effects. The velocity disturbance generated at the 
interface at position x due to N disks, each one centred at 
position x„ and subjected to forces f„ with n = 1,... N, 
is given by v(x) = Y.n=i 2 ^^,+J)||x-x„|| 


2. Narrow gap between two fluid layers 


In this section we study the case where the cell in¬ 
teracts with the substrate only through discrete sites 
(disks) corresponding to the force-generating elements 
while the remaining part of the interface is allowed to 
slip (v(°^ ^ v*^®^) without friction. To this end we con¬ 
sider a narrow gap, a quasi-2D film of height h and small 
viscosity r]a, separating the two fluids (see Fig®. For 
simplicity we shall treat the film as a 2D fluid |(^. 1^. 1^ . 
The coupled Stokes equations now read 

= 0; with V • =0; z > 0 

? 7 aVyV^“) - -I- = -f(5(x); 


with VII 


r(a) _ 


= 0; z = 0 


Ty.V^v^®) - Vp(®) = 0; with V • v^®) =0; z < 0 (A3) 

where indicate velocity and pressure for the in¬ 

tracellular flow; and v*^®\ velocity and pressure for 
the substrate flow as before. p^°''> are the 2D veloc¬ 
ity and pressure in the film, so rja is a two-dimensional 
viscosity and the ratio tja/{qc + Vs) has the dimension of 
a length [sTI. IssMssjl . The terms and —TygCl^vi®) 

describe the stress exerted by the cell and the substrate 
on the film. Equal and opposite stresses are exerted by 
the film on the upper and lower bulk fluids, which en¬ 
ter as boundary conditions for the first and the last line 

of Eq.(|A3|. 

We impose the continuity of the velocities at the loca¬ 
tion of the disks, v(‘’i|disk = v(®)|disk = v(“i|disk- This 
way, using the ansatz vl®)(q) = S(q)e^‘^; v('’’)(q) = 
C(q)e“^'^ ; (q) = A(q) we get the equation in Fourier 

space for the film. Taking the inverse Fourier transform, 

• f. The integral 


= f q 


[l-q(8)q] 


^ ('') I (2-7r)^^ [q’^rja+qi-nc+rje)] 

can be calculated analytically, the result is expressed in 
terms of special functions. The limit of negligible film 
viscosity can be obtained directly by setting ? 7 a = 0 in 
the above expression. The flow = v^°''^x resulting 
from a force in the same direction f = fx is again of the 
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form (r) 


{ria+Vc)r 


/■ 



at the interface at the position of a disk as (for simplicity 
we remove superscripts) 


V||(x,t) 


(Pqdn e-fo're-*f2t 
(27r)2 q[rj^ + 2i 


[i5(n + w) — (5(fl 


a;)]. 


Performing the integrals, we can recast the resulting ex¬ 
pression using Green’s function as 


FIG. 4. (Color online) (a) scheme of a viscous cell of viscos¬ 
ity rjc, lying of a viscous substrate of viscosity qs separated 
by a small gap h of viscosity rja- The disks represent force¬ 
generating elements which are in contact with both the cell 
and the substrate, (b) When subjected to a force, each disk 
drags its surrounding molecules and generates a local flow. 
This flow propagates and affects also other disks far apart. 


^(x,t) = ■ F + (9t[G(")(0 ■ F]; 


(A5) 




(r) = 


27rr 






27rr 


where Sij is the Kronecker symbol. We note that letting 
fj,c ^ 0 one recovers the expression for the semi-infinite 
viscous substrate, and similarly letting ijs ^ 0 one ob¬ 
tains the equivalent expression for the elastic medium. 


3. Equations describing elastic cells lying on 
viscous substrates 

Here we generalise our description to the case of an 
elastic cell characterised by bulk elasticity ^c- The disks 
are arranged with the same geometry considered in Fig HI 
they oscillate around their equilibrium positions of FiglT] 
but are now permanently bound to the cell elastic body. 
The dynamics at the interface (z = 0) is derived as in Sec¬ 
tion ixn In absence of motions along z-axis, the force 
balance is described by = —F(r,z, t) 

where represents the deformation field in the elas¬ 
tic medium. The decomposition of into components 
that are parallel, and orthogonal, ui‘^\ to the in¬ 
terface is defined as in Section ixn The velocities in the 
upper medium (cell) and in the lower medium (substrate) 
coincide at the locations of the disks 

Idisk = Idisk- (A4) 

Once the force-distribution is known, these equations de¬ 
termine both deformation in the cell and the flow gener¬ 
ated in the substrate at the location of the disks. 

Once again to proceed we define Fourier trans¬ 
forms, now with respect to both space and time, 
5 (x,t) = /iga^e-*'i‘-e-*f'‘g(q,0) for any function 

g. From Ea. (|A4p . —To ob¬ 

tain the z-dependence of the solutions we proceed as 
in Section IATI posing V||®^ (q, z, O) = S(q, 0)e'J^ and 

U||'^^ (q, z, O) = C(q, 0)6“'^’’’. Therefore at the disk lo¬ 
cations —iOC(q, O) = S(q, O), resulting in S(q, O) = 
F(q,0)/[q(77s -biff)]. 

As we are interested in the oscillatory behaviour of 
the cell, we consider a point-like force which oscillates 
with frequency w, F(x, t) = fo(5(x) sin(a;t). As a result 
F(q, fl) = fo||[i5(H -b w) — (5(H — w)]. Using this we can 
take the inverse transform and obtain the velocity held 


Appendix B: Dynamics of the four disks system 

We focus on the system made up of two oscillating 
units, labeled with a = A, B, each unit consisting of a 
single dipole discussed in the text, see FiglTJ Although 
in the main text we restricted to a one-dimensional force 
distribution, in the following section we derive the equa¬ 
tions in vectorial form valid for generic force distribu¬ 
tions. 


1. Derivation of the equations describing the 
crawler’s dynamics 

The dynamics at the interface of two fluids with neg¬ 
ligible inertia and viscosities gc and gs is instructive to 
highlight similarities and differences between the locomo¬ 
tion of a soft object and the swimming in bulk viscous 
huids at low Reynolds number [s^. 

We begin with noting that the motion of a disk lying at 
the interface of the two fluids can be very complicated : 
a disk straddling the interface can cause deformations or 
instabilities of the interface which in turn affect the disk 
motion. A simple approximation consists in neglecting 
this aspect of the interface dynamics. Such an approxi¬ 
mation is valid if the surface tension is constant and suf¬ 
ficiently large to prevent any interface deformation. This 
simplified version of the problem has been studied and 
solved by Ranger [^. The analogous problem for the 
sphere was discussed more recently by Pozrikidis (s^ . 

A soft cell migrating over a substrate has peculiar 
mechanical features, which we now explain. The cell- 
substrate interface divides the system into two, individu¬ 
ating intracellular and extracellular forces and flows. In¬ 
tracellular forces originate from the active behaviour of 
the cell and the cytoskeleton. They are responsible for 
cell oscillations but do not contribute directly to the cell 
migration, due to the confining effect of the cell bound¬ 
ary. In our model, only substrate-mediated interactions 



















are responsible for net translational motion of the cell rel¬ 
ative to the substate, which occurs thanks to the stress 
transmitted via the substrate. However, both the dis¬ 
sipation associated to the intracellular forces and that 
associated to the extracellular forces must be included as 
local contributions. This is at the origin of the biphasic 
behaviour as a function of the substrate viscosity in our 
model. 

To see this, we consider the four-disks system of FiglT] 
We write the force balance for the disks as 


F 


a 

n 


C(x“ - v(x“)) -h A° 


(Bl) 


where the new terms A“ are Lagrange multipliers that 
prevent net motion of the centre of the collection of disks 
with respect to the cell frame, representing e.g. the in¬ 
teraction with the cell boundary which is not explicitly 
included in our model. As in the main text, C = Cs + Cc 
is the total drag coefficient ; x“ is the velocity of parti¬ 
cle n ; v(x“) is the flow at particle {a,n) due to all the 
remaining force centres. 

We identify the non-local, intracellular-mediated, force 
acting on particle (Q:,n) : C(x“) := ^cv(x“). This force 
is felt by the disk (a,n) at position x“, due to all the 
remaining particles and represents a non-local effect me¬ 
diated by the intracellular environment. Similarly, we 
identify the non-local, substrate-mediated, force acting 
on particle (a,n) as «S(x“) := ^sv(x“). The meaning 
is analogous: •S(x“) is the force felt by particle (a,n) 
(at position x“) due to all the remaining particles rep¬ 
resenting a non-local effect mediated by the substrate 
(extracellular environment). 

We write Ea. dBlI) inserting these definitions of non¬ 
local intracellular and extra cellular forces and rearrang¬ 
ing terms, Cx“ = F“-|-C(x“)-t-tS(x“) —A“. The Lagrange 
multipliers A“ are determined by requiring : i) no net rel¬ 
ative motion between the centre of the system of oscillat¬ 
ing units and the cell boundary ; ii) that the oscillatory 
dynamics of each dipole da remains unaffected. 

Condition i) is implemented by taking the sum over all 
the particles in the previous equation. By definition the 
system is force-free, so the sum of all the active forces 
driving the particles vanishes, B^n = 0- 

Again by definition, Sa=A B^n = 4c. Condition 

i) reads ELi Ea=A.B C(x“) = Y.l=iY.a=A.B K- To¬ 
gether these result in 


c 


1 c. 

4 (Cs -b Cc) 


E E v(x;) 


(B2) 


In addition, condition ii) implies A“ = A 2 for a = A, B 
so each oscillating unit satisfies 

Cdc = r“-F“ + C[v(x“)-v(x“)]. (B3) 


Analogous expressions can be obtained for the eight disks 
system, see [i^ . 


2. Elastic cells on viscous substrates 

We return to the case of disks bound to an elastic 
cell body. The equation describing the dynamics of one 
isolated disk, centred at position at the interface 

between elastic cell and viscous substrate is given by 
—CsX — Cc[x — -|- f = 0. Here Cs describes the 

substrate viscous drag as before and Cc the coefficient 
relatiM deformation and force applied to the centre of a 
disk lying at the interface (for a disk of radius a, Cc = 
^afj,c). This expression can be readily generalised to the 
case of many disks, labeled as before with indices a and 
n, as -Cs[x“-v(x“)]-Cc[x“-xl'='^^’“-u(x“)]-kF“ = 0. 
v(x“) and u(x“) represent the velocity field in the sub¬ 
strate and the deformation field in the cell at the location 
of disk a,n. From these equations the dynamics of each 
oscillating units a = A, B is obtained as 

do = [F?-F2] -Ed„+v(x?)-v(xg) + |^[u(x?)-u(x-)] 

Ss Ss Ss 

(B4) 

Since the disks are rigidly connected to the elastic 
body, the terms A“ used above to enforce the average 
disk locations with respect to the cell boundaries are not 
needed. The migration speed is obtained by studying the 
motion of tracers lying in the substrate, below the disks, 
at z —>■ 0“. We need four tracers for the four disks sys¬ 
tem and eight tracers for the eight disks. The tracers are 
convected by the flow underneath the disks so the tracer 
dynamics is described by — Cs[x“ — v(x“)] = 0. Thus the 
centre of the collection of tracers moves according to 

a—A,B n—1 

where v(x“) is given by Ea. dASI) . This expression de¬ 
scribes the propulsive speed of the cell. 


3. Force-dipole and quadrupole 


From now on, we restrict ourselves to the one¬ 
dimensional force distribution discussed in the main text. 
Just like in the main text, the disks are centred at 
positions and subjected to forces F“ which satisfy 
/“ := = —F 2 , for a = A, B. The coordinates can be 

written as 


= c- --b(-)"+^ —■ 
2 ^ 2 ’ 


(B6) 


where r = cb — ca is the separation between the centres 
of the two oscillating units, see Fig HI and (—)" = (—1)". 
The quantities c and La are the crawler’s centre and the 
(time-dependent) length of dipole a. With no loss of 
generality we pose La = la + da where la is a constant 
and da is time-varying as in Eq.(|31). 
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The force-multipoles are defined as moments of 
the force distribution. To illustrate this we con¬ 
sider Ea. dBBI) . The moment of order k is := 

Ea=A.BE„=i, 2 -P'“« - c)''- We pose V ■= for 

the dipole, and Q := for the quadrupole, and find 

that they are given by the matrix relation 

(a) = {X Z) ■ ijt) 

showing that the total dipole V of the system is the sum 
of the two individual dipoles. In turn, the two forces 
fA and /b can be expressed as functions of dipole and 
quadrupole moments as 

\ Lb tLb / 

This relation implies that higher moments with 

n > 2, can be expressed as combinations of the two mo¬ 
ments V and Q. 

4. Migration speed 

To begin with, we note that the flow generated at posi¬ 
tion a;“ (the centre coordinate of the generic disk (n,a)) 
by all the remaining disks can be written as 

V{XZ) = Y. Hncm^F^ + E E HncmXi (B9) 

m^n ^^oim—1,2 

where we have introduced Hnamp '■= = 

---g- as obtained in Section ixn 

2'7r(ris+'nc)\x^—x^\ 


where we defined H' := -§?H{r) and H'" := ^H{r). 

c. Instantaneous migration speed and dipole/quadrupole 
moments 

We now briefly discuss how the migration speed is re¬ 
lated to the dipole and quadrupole term of the force- 
distribution. Inserting the expression of /a and fs as 
functions of V and Q obtained from the second relation 
of Eq. (IB8I) we obtain 

^ Ri [iX - XX + 2{L\ + X)—] }• 

(BI2) 

Ea. (IB12p shows that, to leading order in the separation 
r, the sign of the instantaneous migration speed is de¬ 
termined by the quadrupole, consistently with what re¬ 
ported in |26l |. Since here the terms fa and da with 
a = A, B are oscillating quantities with zero mean, the 
same holds for dipole and quadrupole V, Q. As a con¬ 
sequence, also the first term of Ea. (IB12p oscillates with 
zero mean and does not contribute to the net migration 
speed. A net contribution comes instead from the re¬ 
maining terms in Eq. (IB12I) which depend again on dipole 
or quadrupole but are more involved. Ref. [1^ does not 
report how the average migration speed depends on the 
dipole or quadrupole so a direct comparison with exper¬ 
iments is not yet available in this case. 

d. Approximated expression for the average migration speed 



a. Instantaneous migration speed 


From the force balance, using Ea. (IB2l) and the expres¬ 
sion in Eg. dBQI) we obtain the instantaneous migration 
speed for the system depicted in FiglT]as 


c = ^-[{HiAlB — HiA 2 B + H2AIB — H 2 A 2 B)fB 


+ {HiBlA — Hib2A + H2BIA — ^^2A2b)/a]- (BIO) 


b. Approximate expression for instantaneous migration 
speed 

An approximate analytical expression for the migra¬ 
tion speed is obtained by performing the analogue of a 
multipole expansion, valid for 2r ^ |iA + Lb\, up to 
third order in the separation r 

c^^-^[2H'[LBfB-LAfA] (BII) 

+ A- 'iL\LB)fB — {La + 3TaA3)/a] | 


To obtain the average migration speed we insert the 
parametrisation of Eq. ([3]) in Eq. (IB11|) and take the aver¬ 
age over the period T := ^. In doing so, for example, we 
find that the term Ta/a~ Ag/s contains a combinations 
of and and therefore vanishes when we average 
over the period T as discussed above. For the same rea¬ 
son, the term L^Ja — L^Jb does not contribute. Instead, 
the remaining term gives a finite contribution to the av¬ 
erage as LALB{LAfB - LbJa) ~ [2P -f dAdB]{dAfB - 
dBfA) where ^ neglects terms that average to zero. In 
particular, we find dAfB — dBfA ~ —^CloRaRb sinijj 
where fj := c^b - 4>a- Similarly dAdsldA/s - ^b/a] ^ 
— ^R\R^siii2ip. Hence, to leading order in the multi¬ 
pole expansion the average migration speed in a period T 
is the expression reported in the main text, Eq.(|4|), where 




(B13) 


The last equality holds by virtue of the relation among 
deformation and force ga = {CRa^)/‘^ valid at lowest 
order in ajLa- 












10 


To produce a plot we arbitrarily choose a phase shift 
Ip = 7r/2 between the oscillators. For a disk (p = 
[16a(?7s + ? 7 c )]/ 3 . Moreover, we pose Ra = Rb = R 
and QA = 93 = 9 - So, in the case where oscillations 
are driven by providing fixed oscillation amplitude R 
we obtain Vc = i— \ 4 —auPR?. Instead, if oscilla- 

\Xls~r'nc)T' 

tions are driven by providing fixed force amplitude 9 we 

To distinguish between 
these two cases we further pose Sg := and 

S/j := ^auiPR^ and measure the migration speed respec¬ 
tively in units of ^ or in units of keeping r]c con¬ 
stant and varying r]s see Fig.(2)(a). A different choice 
for the phase shift, tp ^ 0,7r, will introduce corrections 
associated to the term cosip in S. As the amplitudes of 
oscillation must satisfy R < I, these corrections do not af¬ 
fect the qualitative trend of the migration speed although 
they can shift the value of substrate viscosity at which 
the migration is optimal, as well as the maximum speed. 


e. Approximated expression for the average migration speed 
for elastic cells on viscous substrates 

In this case Ea. dBSII depends on v(x) given by Ea. dASI) 
and contains the sum of two terms. The first term re¬ 
sembles what was found for the viscous cell lying on the 
viscous substrate. The second term describes the contri¬ 
bution of the elastic interactions. It is important to note 
that such a term contains an exact time-derivative, of 
periodic functions of period T = ^. Thanks to that, 

the average over one period gives dtdt[G^‘^'> ■ F] = 
Qie) . • F|t=o = 0- Thus for calculating 

the average speed we can approximate v(x) ~ • F 

in Ea. dASI) . The multipole expansion for the elas¬ 
tic cell is identical to the one obtained for the viscous 
cell, Ea. dBll|) (and related Ea. dB12|) l. The difference is 
that for the elastic cell forces fa and deformations da sat¬ 
isfy a relation which follows from Ea. dBdl) . Neglecting, as 
before, the effect of the interactions at this level, Ea. dBdl) 
yields a local relation between forces and displacements 
of the oscillators da = —hda + ^fa where T := — is the 
timescale obtained combining the elastic modulus of the 
cell and the viscosity of the substrate. Note that has 
the dimension of a viscosity, so ujT = is equivalent to 
a viscosity ratio. We obtain different force-deformation 
relations depending on the chosen prescription. In par¬ 
ticular we note the parametrisation given in Eq.Q needs 
to be modified to include the additional timescale T, see 
below. 

Assigning forces. If we assign forces as fa = 
—gaSin^Lot -I- (pa) then displacements are obtained solv¬ 
ing the differential equation for da, whose solution at 
the steady state is da{t) ^ [sin(a;t + (pa) — 

ojT cos(a;t -I- (pa)]- As for a viscous cell, the dipole 
term does not contribute to the mean speed while the 
quadrupole does. Posing again ga = 9 and fa = f for 


a = A,B the result is 

3 gsg^T^l^ cos Ip] sin iP 

~ „2 + 4 ATTCsT^l+uj^r^)^ 

IS ' 

(BI4) 

Note that the term (^) has the dimension of length, 
and describes the oscillation amplitude da- Therefore 
it must satisfy (^) < I- Studying the behaviour as a 
function of gg we recover the non-monotonic trend found 
for the case of fluid cell, peaked around Vs ^ fj- It 
is interesting to note that this migration speed is also 
a non-monotonic function of uj, peaked around w ^ ^ 
(weakly sensitive to the values of ip and to the ratio 
(|^)). This result suggests a cell may adapt its oscil¬ 
lation frequency to the environment in order to optimise 
its speed. To produce a plot as a function of the fre¬ 
quency, in this case we consider the expression of Vc given 
by Ea. (IB14l) for the value ip = tt/ 2 as before, obtain¬ 
ing Vc = Figlllb) of the main text 

shows Vc in unit of the velocity measuring uj 

in units of 

The physical reason beyond the existence of an optimal 
frequency in the case of prescribed forces can be under¬ 
stood using a simple argument . The propulsion (and 
pumping of fluid) are proportional to terms that con¬ 
tain product of forces, fa, and deformations, da, of the 
oscillators. Imposing the force, in presence of viscoelas¬ 
ticity, the deformations take place on typical timescale 
7”. Then: (i) for sufficiently high frequencies, T <C T, 
the oscillator cannot reach the maximum amplitude of 
oscillations. In this regime the allowed amplitude of de¬ 
formation increases at increasing T (i.e. at decreasing 
w); (ii) on the contrary for sufficiently small frequencies, 
T 3> T, the oscillator can reach the maximum ampli¬ 
tude available. So in this regime the amplitude of de¬ 
formation does not vary at increasing T. To calculate 
the average speed we must divide da by T. The result¬ 
ing function in regime (i) can increase at increasing T 
while in regime (ii) only decreases at increasing T (be¬ 
having as ^ constant/r). The optimal frequency lies at 
the crossover between these two regions. The situation 
is different when deformations are prescribed functions. 

Assigning deformations. If we assign deformations as 
da{t) = Rcos{ujt -I- (pa) The forces are also prescribed 
functions of time, directly related to the deformations via 
fa = ^[da + ^] which yields fa{t) = wTsin(wf-|- 

(pa) + cos{ujt (pa)]- Also here only the quadrupole con¬ 
tributes to the average as 

Is W 

Interestingly, also this expression has a maximum as a 
function of the substrate viscosity but only monotonic 
trend as a function of uj. Figl^Kb) of the main text shows 
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the velocity as a function of the frequency for fixed de¬ 
formation. In this case we consider the expression of Vc 


given by Eg. dBlSp for the value ip = 7r/2 as before, ob- 

-1 3 j ^'21'2 

taining Vc = ^ We plot Vc in units of the 

velocity measuring lo in units of T“^. 
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